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This  paper  addresses  issues  related  to  output  feedback  control,  including  sensor  place¬ 
ment,  for  a  model  of  an  air-breathing  hypersonic  vehicle.  The  model  presents  a  number 
of  control  challenges,  in  particular  because  of  strong  couplings  between  the  propulsive  and 
aerodynamic  forces.  Because  of  the  vehicle’s  low  weight,  slenderness,  and  length,  the  ve¬ 
hicle’s  flexibility  has  a  large  impact  on  stability  and  control  of  the  vehicle.  Two  output 
feedback  control  methods  are  developed.  One  applies  reconstruction  of  the  flexible  body 
system  states,  toward  applications  of  state  feedback  control.  The  other  uses  a  robust  de¬ 
sign  that  does  not  rely  on  an  observer  to  ensure  stabilization  and  performance  throughout 
a  given  flight  envelope.  A  rate  gyroscope  and  an  accelerometer  have  been  modeled,  in¬ 
corporating  the  flexible  effects,  and  strategies  for  sensor  placement  have  been  developed 
for  the  hypersonic  vehicle  model  to  enhance  observability  or  to  preserve  certain  system 
structures  that  are  favorable  for  robust  control  design.  Simulation  results  are  provided  to 
demonstrate  the  sensor  placement  strategies  and  output  feedback  control  performances. 

I.  Introduction 

Hypersonic  air-breathing  vehicles  offer  a  promising  technology  for  cost  efficient  and  time  reduced  flights 
in  both  commercial  and  military  applications.  Quick  response  and  space  access  meet  the  Air  Force’s  need 
for  global  strike  capabilities  while  providing  significant  advantages  over  expendable  rockets.  Recent  success 
in  NASA’s  X-43A’s  hypersonic  vehicle  affirms  the  scramjet  engine  technology  for  these  applications. 

Control  designs  are  challenging  because  of  unique  characteristics  of  the  hypersonic  vehicle,  particularly 
the  strong  coupling  between  propulsive  and  aerodynamic  forces.  The  engine  thrust  affects  the  pitching 
moment  due  to  the  underslung  location  of  the  scramjet  engine.  The  length,  slender  geometry,  and  relative 
light  weight  of  the  vehicle  cause  the  vibrational  modes  to  significantly  affect  the  aerodynamic  forces.1,2  In 
addition,  the  model  has  critically  stable  internal  dynamics  and  is  statically  unstable.  A  variety  of  approaches 
have  been  utilized  for  the  general  control  problem  in  this  setting;  for  example,  designing  a  nonlinear  controller 
that  neglects  the  flexibility  effects,3,4  or  using  various  linear  state  feedback  controllers  .5~7  Related  work 
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Figure  1.  Geometry  of  the  hypersonic  vehicle  model  from  Bolender  and  Doman.1 


in  the  field  considers  longitudinal  control  through  linearized  models  and  control  designs,2,8  guidance,9  and 
nonlinear  methods,10, 11  many  with  dependence  on  access  to  full  state  information. 

This  paper  utilizes  linear  control  design  by  considering  two  output  feedback  control  methods.  In  the 
first  approach,  partial  state  information  is  provided  through  the  outputs  and  used  to  reconstruct  full  state 
information,  which  is  then  used  with  a  state-feedback  controller.  The  placement  of  the  sensors  strongly 
affects  the  system’s  observability,  so  sensor  placement  strategies  are  developed  to  increase  observability  or 
to  preserve  certain  pole-zero  structures  that  are  favorable  for  robust  control  design.  Our  results  for  the 
hypersonic  vehicle  follow  traditional  sensor  placement  work  for  flight  vehicles  based  on  observability.  In 
particular,  Van  der  Velde  and  Carignan’s  free-free  flexible  beam  model12  closely  agrees  with  our  sensor 
placement  results  for  the  observer  design.  We  also  follow  Al-Shehabi  and  Newman’s13  sensor  placement 
strategies  and  optimization  methods. 

In  the  second  approach,  a  robust  output  feedback  strategy  that  does  not  rely  on  state  observation  is 
developed  to  overcome  the  limitations  posed  by  linear  observer-based  controller.  The  methodology  employs 
pre-compensation  of  the  unstable  zero-dynamics,  dynamic  extension,  and  a  robust  servomechanism  design 
based  on  time-scale  separation  methods.  The  resulting  controller,  paired  with  a  judicious  placement  of 
the  sensor,  is  shown  to  yield  superior  tracking  performance  and  improved  robustness  versus  observer-based 
design  for  a  large  envelope  of  operating  conditions. 

The  paper  is  organized  as  follows:  Section  II  describes  the  hypersonic  vehicle  model  and  sensor  modeling. 
In  Section  III  the  observer  based  output  feedback  controller  is  developed  along  with  relevant  sensor  placement 
strategies.  Section  IV  discusses  the  development  of  the  robust  output  feedback  design,  including  appropriate 
sensor  placement.  Simulation  results  are  provided  in  Section  V,  and  conclusions  are  given  in  Section  VI. 

II.  Hypersonic  Vehicle  Model 

The  hypersonic  vehicle  model  considered  in  this  paper  is  the  assumed-modes  version14  of  the  model 
originally  developed  by  Bolender  and  Doman.1,15  In  the  cited  references,  the  stability-axis  equations  of 
motion  for  the  longitudinal  dynamics  of  the  hypersonic  vehicle  are  derived  through  Lagrange’s  equations  as 

=  — (T  cos  a  —  D)  —  g  sin(d  —  a) 

=  — —  (— T  sin  a  —  L)  +  Q  +  cos  6  —  a 
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Figure  2.  Mass  normalized  flexible  mode  shapes  as  the  fuel  level  is  varied  from  empty  to  50%  90%  -,  and 

100%  full. 


where  to  is  the  vehicle  mass,  Iyy  the  moment  of  inertia,  g  the  acceleration  of  gravity,  and  Q  and  uy  are  the 
damping  factor  and  the  natural  frequency  of  the  z-th  flexible  mode,  respectively.  The  aerodynamic  forces 
enter  (1)  as  lift  L,  drag  D ,  thrust  T,  pitching  moment  M  about  the  body  y-axis,  and  the  generalized  elastic 
forces  Ni,  N2,  IV3.  The  control  inputs  up  =  [<5e,  Sc,  0,  Ad\T  appear  implicitly  in  (1)  through  the  aerodynamic 
forces,  where  5e  is  the  elevator  deflection  angle,  Sc  is  the  canard  deflection  angle,  0  is  the  fuel  to  air  ratio, 
and  Ad  is  the  diffuser  area  ratio.  The  plant  model  considered  in  this  study  comprises  the  eleven-dimensional 
state  vector  xp  =  [Vt.  a,  Q,  h,  9, 771, 771 , 772, 772, 773 , 773]^.  The  first  five  states  are  referred  to  as  xr,  the  rigid  body 
states:  vehicle  velocity  Vt,  angle  of  attack  a,  pitch  rate  Q,  altitude  h,  and  pitch  9.  The  last  six  states  are 
referred  to  as  x/,  the  flexible  states,  which  consist  of  three  vibrational  displacements  771,  772,  773  and  their 
time  derivatives.  A  sketch  of  the  vehicle  geometry  is  given  in  Figure  1. 

The  assumed  modes  version  of  the  hypersonic  vehicle  model14  considers  the  vehicle  as  a  single  flexible 
structure  instead  of  clamped  beams.  This  model  removes  the  coupling  between  rigid  body  and  flexible  body 
dynamics,  so  that  the  flexible  modes  are  orthogonal  to  the  rigid  body  modes.  The  flexible  effect  interaction 
occurs  only  through  the  aerodynamic  forces  in  (1).  Only  the  first  three  modes  are  considered  in  this  study. 
The  vibrational  model  is  defined  as  a  traditional  free-free  beam16  with  mass  normalized  mode  shapes 
shown  in  Figure  2.  The  damping  coefficient  Q  =  0.02  is  constant,  while  the  modal  frequencies  ay  change 
with  the  operating  conditions.  The  most  noticeable  effect  influencing  mode  shapes  and  frequencies  is  the 
mass  of  the  vehicle,  which  decreases  as  fuel  is  consumed.  A  fuel  consumption  model  is  not  yet  developed 
for  the  aircraft,  so  the  fuel  mass  at  trim  condition  is  held  constant  throughout  simulation.  It  is  important 
to  observe  that  while  the  shapes  do  not  change  much  with  fuel  consumption,  there  are  significant  variations 
in  node  locations  and  modal  frequencies.  The  four  sets  of  curves  in  Figure  2  show  the  mode  shape  shifting 
behavior  as  fuel  is  consumed.  A  sample  of  change  in  modal  frequency  and  node,  anti-node  and  center  of 
gravity  locations  for  various  fuel  conditions  is  listed  in  Table  1. 

Table  1.  Selected  first  mode-shape  frequencies  and  node/ ant  i-node  locations  as  vehicle  mass  changes. 


Available  Fuel 

0% 

50% 

90% 

100% 

TO  (slug) 

9356.58 

14786.39 

19130.24 

20216.21 

cu  1  (rad/s) 

22.77672 

21.17379 

20.3468 

20.175 

xcg  (ft) 

53.1008 

53.8218 

54.1039 

54.1555 

^noded  (ft) 

27.4 

30.4 

31.9 

32.2 

node,l  (ft) 

50.5 

51.1 

51.6 

51.7 
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A.  Rate  Gyroscope  Sensor  Model 


The  pitch  rate  of  the  aircraft  can  be  measured  by  a  rate  gyroscope  sensor;  however,  the  rigid  body  pitch 
rate  Q  will  be  affected  by  the  angular  velocity  contribution  of  the  flexible  beam.12,1'  The  measured  pitch 
rate  of  the  aircraft  Qfj,  as  sensed  by  the  jth  rate  gyroscope,  is  modeled  as 


Qf,j 


q-E 


d(f)i  (Xrg.j)  ^ 

,  'Hi 
ax 


(2) 


where  xrgj  is  the  location  of  each  j  rate  gyroscope  matched  with  spatial  derivative  for  each  ith  mass  nor¬ 
malized  mode  shape  in  Figure  2. 


B.  Accelerometer  Sensor  Model 


Many  aircraft  designs  use  a  self-aligning  vane  or  flow  direction  sensor  to  measure  the  angle  of  attack;  how¬ 
ever,  these  traditional  sensors  will  not  suffice  for  hypersonic  flight.  Instead,  we  estimate  the  angle  of  attack 
trajectory  from  measurements  of  normal  acceleration.1'  The  body  axis  normal  acceleration  at  the  center  of 
gravity,  unaffected  by  the  flexible  states,  is  denoted  nz.  Because  our  aircraft  model  only  considers  longitu¬ 
dinal  motion,  we  neglect  rolling  motion  and  normalize  the  acceleration  to  g-units.  This  is  implemented  in 
simulation  with  the  stability  axis  equation: 

nz  —D  sin  a  Lcosa  . 

—  = - b  costh  (3) 

g  mg  mg 

Away  from  the  center  of  gravity  at  a  distance  of  x^  from  the  nose,  each  kth  accelerometer  is  affected  by 
pitch  acceleration,  which  contributes  to  the  measurement  of  rigid  body  normal  acceleration,  denoted  nz y. 


^z,fc  '"z  ( %cg  %a,,k)Q- 

9 


The  flexible  states  affect  normal  acceleration  as  follows: 


(4) 


'"•z  ,f,k 


'"z  ,k 


n 

+ 


q 

y  i= i 


(5) 


Currently,  normal  acceleration  is  not  used  as  a  state  within  the  aircraft’s  equations  of  motion,  so  we  add  (5) 
as  an  output  of  the  plant  linearization. 


III.  Observer  Based  Output  Feedback  Control 

The  first  approach  we  pursue  for  output-feedback  control  is  to  use  a  classical  design  based  on  the  sep¬ 
aration  principle.  For  this  case  study,  a  modified  version  of  the  controller  developed  in  Ref.  5  will  be  used 
as  a  baseline  certainty-equivalence  controller,  to  be  paired  with  an  asymptotic  observer  employed  for  state 
reconstruction.  We  begin  by  presenting  a  linearized  plant  model  and  then  move  on  to  the  design  of  the 
controller. 

A.  State  Variable  Representation 

Previous  control  work5,6  for  the  Bolender  and  Doman  hypersonic  vehicle  model  has  developed  the  framework 
for  the  following  discussion  of  the  state  variable  representation.  As  opposed  to  Ref.  5,  however,  the  diffuser- 
area  ratio  Ad  is  not  used  for  control  design,  and  it  is  therefore  kept  constant  to  the  value  A<j  =  1.  Moreover, 
an  additional  control  effector  (canard)  contributes  to  the  vehicle  control  authority.  The  eleven-state,  three- 
input  nonlinear  system  (1)  is  linearized  about  a  trim  condition  using  numerical  tools.  The  resulting  linearized 
model  is  represented  in  the  familiar  form  as 

Xp  —  ApXp  -|-  BpUp 

Dp  —  CpXp  T  DpUp 

Zp  HpXp ,  (6) 
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State 

Input 

Vt  =  7846.4  ft/sec  m  =  0.54309  ft 

a  =  2.9883  deg  rj2  =  —0.10249  ft 

Q  =  0  deg/sec  ??3  =  —0.034697  ft 

h  =  85000  ft  iji  =  0  ft/sec,  i  =  1,  2,  3 

0  =  2.9883  deg 

5e  =  6.5946  deg 

Sc  =  —5.6041  deg 

0  =  0.20739 

Table  2.  Trim  condition  for  the  linearization  of  the  hypersonic  vehicle  model.  The  vehicle  is  at  50%  fuel 
condition,  with  vehicle  mass  m  =  14786  slug  and  center  of  gravity  xCg  =  53.8218. 


where  the  input  vector  reads  as  up  =  [5e,  5C,  (f)]T .  The  measurement  output  is  given  yp  =  [Vj,  ct,  Qf,j,  h,  nz,f,k] 
(note  the  feed-forward  matrix  Dp  resulting  from  the  normal  acceleration  sensor),  while  the  output  to  be 
controlled  is  given  by  zp  =  [Vt,  7],  where  7  :=  9  —  a  is  the  flight  path  angle  (FPA).  The  trim  condition  used 
in  this  study  to  derive  the  model  (6)  is  given  in  Table  2  The  reference  commands  r  =  [Lref57ref]T  to  be 
tracked  in  this  setup  consist  of  step  functions,  smoothed  by  a  second-order  pre-filter  chosen  so  that  a  1000 
ft/s  change  in  velocity  Vt  settles  after  approximately  100  seconds,  and  a  1°  change  in  flight  path  angle  7 
transitions  after  approximately  2  seconds.  The  reference  model  with  state  xm  G  R2  is  written  as 

•Km  T  BmT 

—  HmXm  ,  (7) 

whereas  the  tracking  error  is  defined  as  e  =  Mmxm  —  Hpxp.  The  model  is  then  augmented  by  the  integral 
error  xi  =  f*  e(r)dr.  The  combined  system,  with  states  x  =  [xp,  x^,  xj  ]T  and  inputs  u  =  up,  is  written  as 

x  =  Ax  +  Bu  +  Gr 

y  =  Cx  +  Du 

e  =  Hx,  (8) 

where  the  matrices  A ,  B,  G,  C ,  D ,  and  H  are  defined  appropriately  in  accordance  with  the  description 
above.  Note  that,  as  opposed  to  Ref.  5,  the  information  available  for  feedback  is  the  measured  output  y. 

B.  Output-Feedback  Control  Design 

Considering  control  design  separately  from  state  estimation,  the  feedback  gains  are  calculated  via  Linear 
Quadratic  Regulator  (LQR)  methods.  The  performance  index  J  of  the  system  is  constructed  to  represent 
the  weighted  sum  of  error  and  input  energy  in  the  following  form: 

1  f°° 

J  =  -  /  ( ztQz  +  uT Ru)dt.  (9) 

2  Jo 

The  selection  of  the  weights  has  been  performed  using  the  implicit  model-following  approach  of  Ref.  5.  The 
optimal  feedback  control  law  reads  as 

u  =  —Kxx  +  Krr  (10) 

where  the  feedback  gain  Kx  is  computed  by  solving  the  continuous  steady-state  Riccati  equation  associated 
with  (8)  and  (9),  and  the  feed- forward  gain  Kr  by  imposing  the  equilibrium  condition  in  steady-state. 

The  physically  available  measurements  of  the  aircraft  only  comprise  three  of  the  state  variables,  so  the 
remaining  states  must  be  estimated  through  the  construction  of  an  observer.  Note  that  only  the  states  of 
system  (6)  need  to  be  estimated,  as  the  integral  error  can  be  reconstructed  directly  from  estimates  of  xp. 
Consequently,  the  expression  of  the  dynamic  observer-feedback  controller  is  given  as 

xp  —  Apxp  Bpup  L(yj  Gpxp  Dpup ) 

u  =  —Kxx  +  Krr,  (11) 

where  x  =  [xp,  xjn,  xj ]T .  The  gain  L  is  calculated  by  solving  a  dual  of  the  LQR  control  problem  (LQE). 
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Eigenvalue 

Observability  Condition 
(sensors  without  flexibility) 

Observability  Condition 
(sensors  with  flexibility) 

-0.49020  -  21.4999i 

1.0506  x  10s 

5.6480  x  105 

-0.49020  +  21.49997 

1.0506  x  10s 

5.6480  x  105 

-1.0771  -  53.71457 

6.6829  x  10s 

2.8714  x  106 

-1.0771  +  53.71457 

6.6829  x  10s 

2.8714  x  106 

-2.1826  -  108.69277 

1.2310  x  109 

4.1545  x  106 

-2.1826  +  108.69277 

1.2310  x  109 

4.1545  x  106 

Table  3.  Results  of  PBH  observability  tests  for  the  system  without  rate  gyros  and  accelerometers  and  observ¬ 
ability  test  with  a  gyroscope  placed  at  1  ft  and  an  accelerometer  placed  at  99 ft.  The  values  are  associated  with 
the  flexible  states. 


C.  Sensor  Placement 


As  it  is  evident  from  equations  (2)  and  (5),  measurements  of  pitch  rate  and  normal  acceleration  are  affected 
by  vibrations  generated  by  the  flexible  dynamics.  The  modeling  of  sensor  vibrations  for  rate  gyroscope  and 
accelerometer  represented  a  key  development,  since  using  simple  rigid  body  sensor  model  in  the  hypersonic 
vehicle  state  observer  proved  unsuccessful  in  simulation  for  full  state  reconstruction.  The  observability  of 
the  flexible  states  system  is  affected  by  placement  of  the  sensors.  Since  the  flexible  states  are  used  in  the 
certainty-equivalence  controller  (10),  the  sensors  should  not  be  placed  at  locations  that  render  the  occurrence 
of  the  flexible  states  negligible.  For  this  reason,  the  role  of  sensor  placement  strategy  is  crucial  for  output 
feedback  design.  A  well  placed  sensor  must  contain  vibratory  components  that  can  be  used  for  more  accurate 
state  reconstruction. 

The  Popov-Belevitch-Hautus  (PBH)  test18  is  used  to  determine  system  observability.  Traditional  con¬ 
trollability  and  observability  tests  involving  the  rank  of  the  controllability  and  observability  matrices  proved 
to  be  unsuitable  for  this  study,  as  the  plant  model  is  poorly  conditioned  numerically.  The  PBH  test  is  more 
complete  and  results  in  a  degree  of  observability,  related  to  the  ratio  between  the  smallest  nonzero  singular 
value  and  the  largest  singular  value  (condition  number)  of  a  certain  matrix.  Specifically,  the  test  states  that 
the  system  is  observable  iff  for  all  eigenvalues  A;  of  A, 


rank 


/ 

1 

1 

^7 

A? 

\ 

V 

C 

) 

=  n. 


(12) 


The  test  has  been  implemented  using  numerical  tools  to  verify  the  observability  of  the  system  and  its  dual 
was  used  to  verify  controllability.  In  the  following  sections  we  will  discuss  the  placement  strategies,  which 
are  defined  by  the  performance  index 


J  =  max(cond(PBH(A,  Ap,  Cp))).  (13) 

For  a  certain  sensor  placement,  we  compute  the  eleven  condition  numbers  of  the  PBH  observability  matrices. 
The  performance  index  J  is  set  to  the  worst  condition  number  (observability  of  the  third  flexible  state  773). 
The  larger  the  condition  number,  the  worse  the  degree  of  observability  of  a  given  mode. 


1.  Independent  Rate  Gyroscope  and  Accelerometer  Placement 

Inspection  of  equation  (2)  shows  that  a  rate  gyroscope  sensor  placed  at  a  dominant  mode  shape  anti-node 
will  minimize  the  contribution  of  the  flexible  term,  setting  Qf~Q.  This  is  a  poor  choice  for  an  observer 
because  the  rate-gyroscope  measurement  does  not  contain  sufficient  information  from  the  flexible  states.  This 
location  produces  an  ill-conditioned  observability  matrix  (for  example,  at  xrg:\  =  51.1  ft  in  Figure  3(a)). 
The  system  outputs  will  contain  richer  information  from  the  flexible  effects  if  the  rate  gyroscope  sensors  are 
placed  at  locations  of  high  ^  values,  such  as  at  mode  shape  nodes  for  example,  at  xTg,i  =  30.4  ft  in  Figure 
3(a)).  This  effect  is  seen  in  Table  3,  where  observability  matrices  have  maximum  order  of  106  for  a  well 
placed  sensor  suite.  The  last  column  in  Table  3  represents  the  improved  observability  condition  numbers  for 
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No.  of  Sensors 

J 

Rate  Gyroscopes 

Accelerometers 

2 

4.1373  x  106 

(Trg?i  —  1.0 

xa,i  =  99.0 

3 

2.8903  x  106 

Xrg.l  =  5.4 

xa,i  =  66.5,  xa,2  =  11-9 

3 

4.3852  x  106 

xrg  i  —  2.7^  xrg52  —  97.3 

xa,i  =  F0 

4 

2.6488  x  106 

xrg)i  =  76.8,  xrg:2  =  14.3 

®a,i  =  88.7,  xa>2  =  12.0 

Table  4.  A  sample  of  results  for  simultaneous  placement  of  multiple  rate  gyroscopes  and  accelerometers. 


well-placed  rate  gyroscope  and  accelerometer  sensors.  The  observability  improvement  can  quickly  be  seen 
in  the  last  row,  which  is  initially  of  order  109  and  improves  to  order  106  after  careful  sensor  placement. 
Figure  3(a)  shows  a  single  rate  gyroscope  sensor’s  effect  on  observability,  plotting  J  in  (13)  versus  physical 
location  xrgj.  The  effect  of  placing  a  single  accelerometer  on  observability  of  the  flexible  stated  is  shown 
in  Figure  3(b).  The  peaks  of  higher  unobservability  index  correlate  with  low  ^  values  evaluated  along  the 
mode  shapes. 


Figure  3.  Rate  gyroscope  and  accelerometer  sensor  placement  effects  on  observability  at  50%  fuel  condition. 


2.  Simultaneous  Rate  Gyroscope  and  Accelerometer  Placement 

Placing  multiple  rate  gyroscopes  and  accelerometers  on  the  aircraft  can  best  be  accomplished  through  the 
application  of  nonlinear  constrained  optimization.  Again,  the  performance  index  for  this  optimization  is  the 
maximum  observability  condition  number  of  the  plant  from  (13).  A  quasi-Newton  algorithm  is  used  for  the 
optimization,  consisting  of  Hessian  update  and  line  search  procedures.  The  locations  are  constrained  between 
the  nose  and  tail  of  the  aircraft;  further,  shooting  methods  divide  the  aircraft  into  ranges  for  each  sensor  to 
be  placed.  The  mode  shapes  were  provided  with  resolution  of  0.1  ft,  so  a  branch  and  bound  algorithm  must 
also  be  applied  to  account  for  these  discrete  steps.  A  sample  of  the  performance  index  results  are  listed  in 
Table  4. 

The  sensor  suite  chosen  for  observer  simulations  in  Section  D  is  a  rate  gyroscope  placed  at  1  ft,  and 
an  accelerometer  placed  at  99  ft.  This  choice  was  made  through  the  aforementioned  reasoning  for  sensor 
placement  strategy,  as  well  as  the  numerical  optimization  results  contained  in  Table  4.  Adding  redundant 
sensors  would  reduce  the  effect  of  sensor  failure,  increase  controller  reliability  and  increase  observability;  but 
this  suite  of  two  sensors  will  provide  the  baseline  for  simulations  of  observer  based  output  feedback  control. 

D.  Simulation  Results  of  Observer  Based  Output  Feedback  Control 

The  nonlinear  equations  of  motion  (1)  have  been  implemented  in  Simulink  for  simulation  with  the  observer- 
based  controller  designed  for  the  linearized  plant.  We  consider  a  default  fuel  condition  as  |  full,  with  trim 
condition  given  in  Table  2.  The  flight  path  angle  reference  7ref  is  increased  by  1°  from  the  trim  conditions, 
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(b)  Flight  Path  Angle 


Figure  4.  Reference  trajectories:  The  vehicle  velocity  is  commanded  constant,  and  the  flight  path  angle  is 
commanded  1°  up-and-down. 


then  returned  to  trim  after  2  seconds,  while  the  vehicle  velocity  Vt  is  commanded  to  remain  constant  at  the 
trim  value.  The  results  of  the  simulation  in  nominal  operating  conditions  (when  the  observer  model  matches 


(a)  Pitch  Rate 


(b)  Angle  of  Attack 


Figure  5.  Pitch  rate  (a)  and  angle  of  attack  (b)  showing  output  feedback  control  results  using  a  gyroscope 
placed  at  1  ft  from  the  nose  and  an  accelerometer  placed  at  99 ft. 


the  linearized  vehicle  dynamics)  reveal  a  somewhat  limited  capability  of  the  observer-based  controller  to 
achieve  tracking  of  the  reference  trajectory.  In  particular,  Figure  4,  showing  the  controlled  outputs  alongside 
their  commanded  reference,  reveals  that  stability  is  maintained  in  closed-loop  and  that  the  tracking  error 
converges  to  zero  asymptotically  after  a  sizable  transient.  Figure  5(a)  shows  the  rate  gyroscope  measurement 
Qf ,  the  actual  state  trajectory  Q,  and  the  reconstructed  state  trajectory  Q.  The  angle  of  attack  estimate 
a  is  also  compared  with  its  actual  trajectory  a  in  Figure  5(b).  The  remaining  states  are  successfully 
reconstructed  with  similar  results.  The  linear  observer-based  controller,  however,  shows  poor  robustness  to 
plant  variations.  Commanding  a  typical  change  in  velocity  of  1000  ft/sec  results  in  inaccurate  state  estimates 
that  destabilize  the  system,  as  seen  in  Figure  6.  Such  a  change  invalidates  the  linearization  to  the  point  that  a 
new  linearization  for  the  observer  must  be  computed.  The  closed-loop  system  is  also  sensitive  to  parameter 
variations,  namely  to  changes  in  fuel  conditions.  The  baseline  linear  observer  designed  for  a  \  full  fuel 
condition  is  simulated  with  a  plant  trimmed  at  empty  fuel  condition  in  Figure  7(a).  The  state  estimates  do 
not  converge  to  the  true  state  values  due  to  this  system  uncertainty,  and  the  inaccurate  estimates  destabilize 
the  system.  Instability  also  occurs  for  a  simulation  with  a  full-fuel  model,  as  shown  in  Figure  7(b). 


IV.  Robust  Output  Feedback  Control 

The  poor  robustness  exhibited  by  the  linear  observer-based  controller,  which  could  not  be  mitigated  by 
a  different  selection  of  the  weights  of  the  LQR/LQE  problem,  motivated  a  further  investigations  into  robust 
output  feedback  methods  that  do  not  rely  on  an  observer.  By  writing  the  linearized  system  in  normal  form 
we  observe  that  the  unstable  zero-dynamics  with  respect  to  the  I/O  pair  (<p:Se)/( Vt,  h)  can  be  compensated 
by  means  of  an  auxiliary  robust  dynamic  controller  which,  driven  by  the  output  Qf ,  commands  the  input 
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7600. 


-Vt  -  Output 
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(b)  Flight  Path  Angle 


Figure  6.  Observer  failure:  The  baseline  reference  trajectory  is  modified  so  that  velocity  is  commanded 
to  increase  by  1000  ft/s.  The  flight  path  angle  reference  trajectory  is  unmodified  from  the  command  of  1° 
up-and-down. 


(a)  Empty  Fuel  Tank 


Figure  7.  Observer  failure:  Baseline  observer  designed  for  half  full  fuel  tank  and  (a)  Simulated  for  an  empty 
fuel  tank  plant,  (b)  Simulated  for  a  full  fuel  tank.  Gyroscope  placed  1  ft  from  the  nose  and  accelerometer 
placed  at  99  ft. 


5C.  The  zero-dynamics  include  the  unstable  pitch  dynamics  and  the  stable  dynamics  of  the  flexible  states. 
The  auxiliary  controller  is  designed  to  be  robust  with  respect  to  parameter  variations  due  changes  in  the 
fuel  condition.  Since  the  position  of  the  rate  gyro  sensor  affects  the  position  of  the  transmission  zeros  of  the 
zero  dynamics,  sensor  placement  plays  a  major  role  in  the  implementation  of  the  auxiliary  controller.  After 
the  zero-dynamics  have  been  compensated,  servomechanism  theory  is  employed  to  design  a  partial-state 
feedback  controller  achieving  robust  tracking  of  reference  trajectories  for  the  vehicle  velocity  and  altitude. 
Time-scale  separation  and  composite  control  methods  are  used  to  alleviate  the  complexity  of  the  design. 

A.  Linearized  Vehicle  Model 

The  dynamics  of  the  flexible  vehicle  model  linearized  about  a  trim  condition  are  re-written  as 

xr  =  An{n)xr  +  A12(n)xf  +  Bi(fx)u 
Xf  =  A2l(y,)xr  +  A22(fi)xf  +  B2(y,)u 
y  =  ClXr  +  C2{n)xf  ,  (14) 

where  xr  is  the  rigid-body  state,  Xf  is  the  flexible  state,  u  =  [(p,Se,Sc]T  is  the  control  input,  and  y  = 
[Vt,  h,  Q f]T  is  the  measured  output.  Measurements  of  the  normal  acceleration  have  not  been  included  in  the 
output  vector,  due  to  their  sensitivity  on  the  operating  conditions.  The  vector  p  G  W  denotes  uncertain 
parameters  due  to  the  trim,  the  position  of  the  rate  gyro  sensor,  and  other  plant  parameters  such  as  the 
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fuel  level.  The  actual  value  of  y  is  not  known,  but  y  is  assumed  to  range  over  a  given  compact  set  V  C  W\ 
which  depends  on  the  flight  envelope.  It  is  assumed  that  only  one  rate  gyro  sensor  is  employed.  Finally,  the 
output  to  be  controlled  is  chosen  as  y  =  [Vt,  h\T ,  where  the  altitude  h  replaced  the  flight-path  angle  7.  The 
following  facts  hold  for  the  plant  model,  for  any  y  E  V: 


1.  The  system  (14)  has  vector  relative  degree  r  =  {1,2}  with  respect  to  the  regulated  output  y  and  the 
control  input  u. 

2.  Any  combination  of  two  different  control  inputs  among  the  three  available  in  u  results  in  unstable 
zero  dynamics  with  respect  to  y.  In  particular,  choosing  u  =  (0,  Se),  the  resulting  8-dimensional 
zero-dynamics  with  respect  to  y  possess  3  pairs  of  complex  conjugate  eigenvalues  in  Re  [A]  <  0  (due 
to  the  flexible  effects)  and  a  pair  of  real  eigenvalues  (due  to  the  pitch/pitch  rate  dynamics,)  located 
symmetrically  with  respect  to  the  imaginary  axis. 


1.  Normal  Form  &  Zero-Dynamics 

Denote  with  waux  £  M  and  yaux  £  M  the  auxiliary  input  and  the  auxiliary  output  of  system  (14),  given 
respectively  by  waux  =  5C  and  yaux  =  Qf,  and  change  coordinates  as  (xr,Xf)  — *  ( x,fj ),  where 

x=[Vt,h,h]T,  f)  =  [riT,f)T,6,Q]T 


and  h  replaces  a  as  a  state  of  the  system.  Thus,  we  obtain  the  system 


x  =  An(n)x  +  A12(n)fj  +  Bn(n)u  +  B12(n)uaux 

f}  =  A2i  (y)x  +  A22  (y)fj  +  B2i(y)u  +  B22(y)uaux 

y  =  CuX 

?/aux  =  C22(y)fj.  (15) 


The  triplet  (An,5n,Cii)  has  the 

following  structure 

^  an 

0  12 

«13^ 

^bn 

bi2^ 

An(y)  — 

0 

0 

1 

,  Bn(y)  = 

0 

0 

\«31 

032 

«33/ 

\p31 

^32/ 

where  the  matrix 

II 

CO  i-* 

bi2\ 
^32  j 

C11  = 


(16) 


is  non-singular.  In  particular,  it  is  assumed  that  a  number  6^  >  0  is  known  such  that  bn(y)  >  6°i  f°r 
y  £  V.  Moreover,  since  it  can  be  shown  that  Im  Ai2  C  Imfln,  the  matrix  A\2  has  the  following  structure 


Ai2(y) 


( *  ... 

0  ...  0 
y*  ...  */ 


where  V  denotes  unspecified  (possibly  non-zero)  entries.  Changing  coordinates  again  using 


z  =  n~  B2i(y)B  1(y) 

and  rearranging  variables,  the  system  is  put  in  normal  form 

z  =  F(y)z  +  G(y)u 

aux  +  Jll(p)x 

x  =  A(y)x  +  B(y)u  + J2i(y)z  +  J22(y)uaux 
y  =  C  x 

!/aux  =  H(y)z  +  J\2(y)x  (17) 
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where  B(y)  =  Bn(/i),  C  =  Cn,  and 


^  an 

dl2 

Ol3^ 

(*  - 

(*\ 

A(y)  = 

0 

0 

1 

5  J21  (ff )  — 

0  . 

.  0 

5  ^22(ff)  — 

0 

VDl 

032 

O33/ 

V*  • 

•  */ 

w 

It  is  well  known  that  the  eigenvalues  of  F(y)  are  exactly  the  transmission  zeros  of  the  original  system 

(between  the  input  u  and  the  output  y .)  Due  to  the  structure  of  the  triplet  (A,B,C),  and  the  fact  that 

Im  [J21  J22]  C  Im5,  the  zero  dynamics  of  the  system  are  given  by  the  trajectories  of  the  controlled  system 

i  =  F{n)z  +  G(y)uaux 

2/aux  =  H(fl)z.  (18) 

It  can  be  verified  that  system  (18)  is  both  controllable  and  observable  for  almost  all  y  G  V,  except  for  certain 
singular  positions  of  the  rate  gyro  sensor,  where  observability  of  some  flexible  state  is  lost. 

B.  Robust  Stabilization  of  the  Zero-Dynamics 

Since  the  triplet  (F,  G,  H )  is  controllable  and  observable,  it  is  possible  to  design  a  dynamic  output  feedback 
controller  to  stabilize  the  zero  dynamics  system  (18).  However,  this  must  be  accomplished  in  spite  of 
the  parameter  uncertainty,  and  thus  a  standard  approach  based  on  state  observation  may  lead  to  poor 
robustness.  A  closer  look  at  the  structure  of  (18)  reveals  that  the  system  in  question  has  relative  degree 
one  and  is  weakly  non-minimum  phase,  in  the  sense  that  all  zeros  are  in  Re  [A]  <  0,  with  a  simple  zero  at 
the  origin.  In  particular,  the  pole/zero  structure  of  the  zero-dynamics  for  the  nominal  trim  condition  is 
depicted  in  Fig.  8.  In  particular,  note  the  pole/zero  interlacing  property  due  to  the  flexible  dynamics,  and 
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Figure  8.  Pole/zero  structure  of  the  controlled  zero  dynamics  (18). 


the  presence  of  the  zero  at  the  origin,  due  to  the  fact  that  the  auxiliary  output  yaux  is  a  linear  combination 
of  pitch  rate  and  derivatives  of  the  flexible  modes.  Changing  the  position  of  the  rate  gyro  sensor  will  affect 
the  position  of  the  complex  zeros.  The  important  issue  of  sensor  placement  will  be  addressed  in  the  next 
section. 

Any  static  output  feedback  of  the  form  uaux  =  — fc?/aux,  k  >  0  results  in  an  unstable  closed-loop,  as  the 
right-hand  side  pole  migrates  towards  the  zero  at  the  origin.  However,  the  positive  feedback  interconnection 
with  a  simple  dynamic  compensator  of  the  form 


Ka{s ) 


=  — K, 


s  +  (3 
s  —  7  ’ 


te  >  0 ,  (3  >  0  ,  7>0 


(19) 


results  in  a  stable  closed  loop  system  if  7  and  (3  are  suitably  chosen  smaller  that  the  magnitude  of  the  unstable 
pole,  and  k  is  large  enough.  Moreover,  if  the  structure  of  the  pole/zero  map  of  (18)  is  preserved  when  y 
varies  within  the  compact  set  V,  the  simple  dynamic  controller  provides  robust  stabilization.  Figure  9(a) 
shows  an  example  of  the  root  locus  of  the  controlled  zero-dynamics,  whereas  Figure  9(b)  shows  the  variation 
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(a)  Root  locus  of  the  compensated  zero  dynamics. 


(b)  Root  locus  of  the  compensated  zero  dynamics.  Variations  in  fuel 
level:  (a)  empty  tank,  (b)  half-full  tank,  (c)  full  tank. 


Figure  9.  Root  loci  of  the  compensated  zero  dynamics. 


of  the  root  locus  due  to  changes  in  fuel  tank  from  full  (a)  to  empty  (c).  It  can  be  noted  that  the  same 
controller  stabilizes  all  three  configurations.  The  robust  compensator  (19)  has  a  state-space  realization 

4  =  7^0  +  {P  +  7)?/aux 

Manx  =  -K&)  ~  «2/aux  (20) 

with  state  £o  £  yielding  the  augmented  internal  dynamics 

za  =  Fa(n)za  +  Jf{p)x 
x  =  Aa(p)x  +  Jf(p)za  +  B(fj)u 

y  =  Cx  (21) 

where  za  =  [fo,zT]T,  and  Fa(/P),  Aa{n),  Jf(/i),  and  J%  are  defined  appropriately. 

1.  Sensor  Placement  for  Robust  Compensation  of  the  Zero-Dynamics 

As  mentioned  above,  placement  of  the  rate  gyro  affects  the  position  of  the  complex  zeros  of  system  (18). 
Moving  the  sensor  from  the  front  of  the  vehicle  towards  the  back  causes  the  zeros  to  migrate  upwards  and 
pass  through  the  poles.  Since  the  sensor  is  placed  at  a  fixed  position,  it  is  desirable  to  obtain  a  pole-zero 
structure  that  remains  favorable  for  robust  stabilization  for  any  fuel  level  and  flight  condition  within  the 
flight  envelope.  By  a  favorable  pole-zero  structure  we  mean  that  the  angles  of  departure  of  the  complex 


Figure  10.  Mass  normalized  flexible  mode  shapes  derivatives  for  a  specific  flight  condition  and  vehicle  param¬ 
eters. 


12  of  22 


American  Institute  of  Aeronautics  and  Astronautics 


Root  Locus 


Root  Locus 


(a)  Root  loci  for  the  compensated  zero  dynamics  for 
two  fixed  fx,  the  rate  gyro  sensor  placed  at  either  25ft 
or  30ft  from  the  front. 


(b)  Root  locus  of  the  compensated  zero  dynamics. 
Variations  include  5  different  flight  conditions  from 
Vt  =  7846ft / s-8846ft /s,  h  =  85000ft-95000ft,  and  2 
fuel  conditions,  50%  and  90%. 


Figure  11.  Root  loci  of  the  compensated  zero  dynamics  &  sensor  placement 


roots  of  the  zero  dynamics  with  the  auxiliary  compensator  are  towards  the  left.  As  seen  on  figure  11(a),  if 
the  sensor  is  placed  passed  a  certain  position,  the  angle  of  departure  for  the  third  flexible  mode  derivative 
turns  towards  the  right.  The  exact  position  of  this  point  depends  on  the  mode  shapes  of  the  flexibility  and 
thus  on  the  fuel  level.  Furthermore,  we  seek  to  minimize  the  effect  of  the  flexible  states  on  the  output  Q/, 
thus  providing  a  better  estimate  of  Q  and  reducing  oscillations  in  our  canard  control  commands  due  to  the 
flexible  states.  In  this  respect,  it  would  be  ideal  to  place  the  sensor  at  common  node  for  all  the  mode  shapes 
of  the  flexible  state  derivatives.  Nodes  correspond  to  pole-zero  cancellations  in  the  zero-dynamics  and  thus 
to  unobservability  of  the  respective  flexible  mode  from  Qf.  However,  as  we  see  from  figure  10,  the  modes 
do  not  posses  a  common  node.  Consequently,  we  have  chosen  to  place  the  sensor  such  that  we  minimize 
a  weighted  norm  of  the  mode  shapes  of  the  flexible  state  derivatives  while  maintaining  conservatively  a 
favorable  pole-zero  structure  of  the  zero  dynamics  for  any  fuel  level  and  flight  condition  within  the  flight 
envelope,  as  described  above. 

C.  Filtered  Transformation 

Once  the  internal  dynamics  have  been  rendered  asymptotically  stable,  a  viable  strategy  to  stabilize  the 
overall  system  (21)  is  to  enforce,  via  u,  a  small-gain  interconnection  using  feedback  from  the  partial  state 
x.  Note,  however,  that  the  state  x  is  not  available  for  feedback.  For  this  reason,  a  preliminary  dynamic 
transformation  is  needed  to  obtain  a  system  that  is  in  the  required  form  for  partial-state  feedback  control. 
To  avoid  differentiating  the  second  component  of  the  output,  a  so-called  filtered  transformation  is  employed. 
For  convenience  of  notation,  let  u  =  [iq,  U2]T  =  [</>,  <5e]T,  y  =  [y i,  y2]T  =  [Vt,  h],  and  augment  the 
dynamics  (21)  with  the  scalar  system 

Ci  =  —  A  Ci  +  «2  (22) 

where  A  >  0  is  a  design  parameter.  Next,  applying  the  change  of  coordinates 

631  ,  ^31^12  —  ^32^11  ^ 

x3  =  -  —  xi  H - - - Ci 

on  on 

and  defining  the  state  of  the  “extended  zero  dynamics”  as  ze  =  [zaT,x 3]T,  one  obtains  the  system 

ie  =  F*(y,)ze  +  JfMy  +  JeMZi 
V  =  Aen(y)y  +  Ae12(y,)^  + J^(y)ze  +  B(y)u 

Ci  =  —  ACi+«2-  (23) 
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In  particular,  the  matrix  Af2(y)  has  the  structure 


A\2{y)  _ 


^31^12— ^32^11 

fell 


where  b3lbl2bii632bl1  ^  0,  since  B(y)  in  (16)  is  non-singular.  As  a  result,  the  further  change  of  coordinates 


^31^12  —  b32bu 


-J{2(n)x  2,  x  =  [x1,x2,£i]: 


yields  a  system  in  the  form 


fe  =  Fe(iu)ze  +  J?(y)y 
x  =  Ae(y)x  +  J2(y)ze  +  Be(y)u 

y  =  Cex ,  (24) 

where  Fe(y)  is  Hurwitz  for  any  y  G  V,  and  the  triplet  ( Ae1Be,Ce )  has  the  following  structure: 


A  (y)  a21 


0  -A 


'&11  f>12 
Be(y)  =00 

\  0  1 


1  0  0 
0  1  0 


It  can  be  shown  that  a23(y)  ^  0  for  all  y  G  V  such  that  the  matrix  in  (16)  is  invertible.  Without  loss  of 
generality,  assume  a23(y)  >  O.Note  that  the  upper  subsystem  (the  new  internal  dynamics  with  state  ze)  is 
interconnected  to  the  lower  subsystem  only  via  the  regulated  output  y .)  Furthermore,  the  partial  state  x 
is  now  available  for  feedback.  Note  that  rearranging  the  order  of  the  state  variables  as  (x2lfi1Xi)  has  the 
effect  of  highlighting  a  triangular  structure  exhibited  by  the  system.  The  triangular  structure  suggests  a 
hierarchical  approach  to  the  control  design.  Specifically,  the  inner  chain  of  integrators  (x2l  £i)  corresponding 
to  the  extended  altitude  dynamics  should  be  controlled  via  u2  =  Se  on  a  faster  time  scale  than  the  dynamics 
for  Vt,  which  is  controlled  primarily  via  U\  =  4>.  In  regard  to  this,  note  that  the  zero  dynamics  of  (24)  with 

respect  to  the  input/output  pair  is  given  by  the  controlled  system 

fe  =  Fe(y)  ze  +  JfM  x2 

X2  =  J22(y)  ze  +  a^2(y)  x2  +  de23(y)  fx 

£i  =  —  A£i  +  u2 

y2  =  x2  (25) 

which  has  relative-degree  two  with  respect  to  the  I/O  pair  u2/y2 ,  transmission  zeros  in  Re[s]  <  0  (the 
eigenvalues  of  Fe(y)),  and  partial  state  (x2,fi)  available  for  feedback.  Once  a  robust  output-feedback 
controller  1C2  for  (25)  has  been  designed,  the  task  becomes  that  of  designing  a  controller  A Z\  for  the  relative 
degree  one-  minimum-phase  system  with  input  ui,  output  xi,  as  shown  in  Figure  12. 

D.  Robust  Servomechanism  Approach 

In  this  section,  we  consider  the  problem  of  designing  a  robust  regulator  for  the  system  (24).  The  problem  is 

that  of  letting  y(t )  asymptotically  track  a  given  reference  signals  r(t).  As  before,  we  let  r(t)  be  the  output 

of  a  decoupled  2nd  order  reference  model  driven  by  step  inputs;  consequently,  r(t)  can  be  generated  as  the 
output  of  the  autonomous  system 

iu  =  Sw 


r  =  Qw 


where  w  =  [wf , w2]T  G  M6,  r  =  [ri,r2]T  G  K2, 


/  \  /  \  0  1  0  \ 

s=\  d  s  ’  «“(  o  o  ’  s,=  -<<  <i  •  «-=(l  0  o),  *  =  1,2 

V  2/  V  V2/  \  o  o  o  / 

and  (Ci,u Jn,i)  are  respectively  the  damping  ratio  and  the  natural  frequency  of  the  i-th  reference  model. 
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Figure  12.  Inner  and  outer-loop  control 


1.  Inner-loop  Design 

Consider  again  the  expression  of  (25),  augmented  with  the  reference  model  for  :r2.  Assuming  that  the  roots 
of  the  Hurwitz  polynomial  ps2(s)  =  s2  +  2(ujns  +  uj^l  are  chosen  to  be  different  from  the  eigenvalues  of  Fe(p), 
there  exists  a  unique  parameter  dependent  solution  II2(/i)  of  the  Sylvester  equation 

n2((j)s2  =  n2(M)  +  . 


Then,  the  change  of  coordinates 


e2  =  X2  -  r2  ,  5e  :=  ze  -  TL(p)w 

representing  respectively  the  tracking  error  and  a  set  of  transversal  coordinates  to  the  zero  dynamics  of  the 
above  system  with  respect  to  e2,  yields  the  error  system 


w2  =  S2w2 

!e  =  Fe^)ze  +  rMe2 

e2  =  J22{p)  ze  +  ae22{p)  e2  +  ae2Z{p)[^  -  R2{p)w\ 

ti  =  -A<£i+w2  (27) 

where 

F2 (/r)  =  _e  t  -  \_J22(fi)D2(p)  +  a22(ti)Q2  —  Q2S2j  . 

a23\tl) 

Note  that  in  order  to  enforce  invariance  of  e2  =  0,  the  state  £i  must  be  steered  to  the  steady  state  trajectory 
1,1  (t)  =  R2(n)w(t ),  which  is  not  available  since  fi  is  not  known. 


2.  Two-time  Scale  Internal  Model  Design 


Assume,  without  loss  of  generality,  that  the  pair  (. R2(n),S2 )  is  observable  for  all  n  €  V.  Then,  the  pair  is 
topologically  equivalent  to  the  pair  (Q2,  S2),  where 


S 


(°  1 

0  0 

VO 


0  > 

1 

-2C  Un) 


As  a  matter  of  fact,  it  is  easy  to  see  that  S2  =  02  1(p)S2Q2(p,)  and  R2(fi)  =  Q 2©2(/i),  where  B2(t)  is  the 
observability  matrix  of  (R2(/j),  S2).  It  should  be  emphasized  that  while  the  pairs  (. R2(/i),S2 )  and  ( Q2,S2 ) 
are  equivalent  representation  of  the  same  autonomous  system,  the  latter  is  independent  of  the  uncertain 
parameter  /x.  This  result  is  exploited  for  the  construction  of  an  internal  model  of  the  uncertain  steady  state 
trajectories  of  (27)  compatible  with  the  condition  e2  =  0. 
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The  goal  of  the  inner-loop  control  is  (i)  to  enforce  by  feedback  an  invariant  subspace  for  the  system  (27), 
in  which  the  tracking  error  is  identically  zero;  (ii)  to  render  the  subspace  in  question  globally  attractive;  (iii) 
to  induce  a  time-scale  separation  between  the  5e-dynamics  and  the  dynamics  of  the  tracking  error  e2.  To 
begin  with,  denote  with  p  >  0  an  adjustable  high-gain  parameter,  and  let  £  =  1/p.  Let  (E2,H2)  be  the  pair 
in  controller  form 


(  0 

1 

0  > 

(o\ 

E2  = 

0 

0 

1 

5  -2  — 

0 

\-0o 

-01 

-02/ 

W 

(28) 


where  the  coefficients  of  the  Hurwitz  polynomial  ps2  (s)  =  s3  +  cr2s2  +  oqs  +  cr0  are  chosen  arbitrarily,  in  such 
a  way  that  the  roots  of  p-^2(s)  are  distinct  from  those  of  the  polynomial  ps(s).  Then,  for  any  p  >  0  there 
exists  a  unique  solution  M2  (p)  of  the  Sylvester  equation 


M2(p)S  —  pE2M2(p)  +  zl2Q2  (29) 

which  is  also  nonsingular.  Post-multiplying  both  sides  of  equation  (29)  by  M2“1(p),  one  obtains 

M2(p)SM2  \p)  =  pE2  +  ^2T2(p) 

where  T2(p)  :=  Q2M2l(p).  It  is  readily  seen  that  T2(p)  is  the  unique  feedback  matrix  that  assigns  to  the 
pair  (pE2,  H2)  the  eigenvalues  of  S,  so  the  expression  of  \k2(p)  can  be  easily  computed  as 


d,2(p)  =  p^2^) ,  d'2(e)  :=  ^-cr0  — ay  +  e2uj2  -a2  +  2^C^n)  •  (30) 

In  particular,  note  that  T2(e)  =  0(1).  A  candidate  internal  model  for  the  steady-state  trajectory  ££(i)  = 
R2(p)w(t)  =  Q2^>2(p)w(t)  is  constructed  as  the  one-parameter  family  of  dynamical  systems 

C2  =  pS2£2  +  “2^1 

v  =  T2(p)£2.  (31) 

Augmenting  the  system  (27)  with  the  internal  model  (31),  changing  coordinates  as 

Ci  :=  Ci  ~  ^2(p)C2 

C2  :=  C2  -  M2(p)02  1(p)w  -  _e  .  52e2  , 

a22.\kJ) 

and  rearranging  terms  after  some  manipulations,  one  obtains  the  system 
5e  =  Fe(p)  ze  +  J{2(p)  e2 

e2  =  J|2(p)5e+ [a|2(p)  +  T2(p)S2]e2  +  a|3(p)Ci+a|3(p)T2(p)^2 
Ci  =  —  ACi  —  d,2(p)  [pS2^2  +  H2Ci]  +  u2 

C2  =  pE2^2-  ^^E2J22(p)ze  +  ^e^[pE2-a^2l]E2e2.  (32) 

Note  that,  for  reasons  that  will  become  clear  in  a  moment,  we  have  kept  the  original  coordinates  Ci,  C2  in 
the  expression  of  the  dynamics  of  Ci-  Since  the  state  variables  of  the  dynamic  extension  and  the  internal 
model  are  available  for  feedback,  and  so  is  the  partial  state  (e2,Ci),  the  control  u2  can  be  selected  as  the 
“feedforward  plus  partial-state  feedback” 

U2  =  [Ai  +  T2(p)S2]er  +  p*2(p)E2£2  -  pk2e2  -  pfcfi  (33) 

where  k  >  0  is  a  scalar  gain  parameter.  Letting  y2  =  col  (e2,  Ci,  C2),  and  dividing  the  corresponding 

dynamics  with  p,  system  (32)  in  closed  loop  with  (33)  is  conveniently  written  as  a  singular  perturbation 

model  in  standard  form 


5e  =  Fe(p)ze  +  Jt2(p)C2X2 

£X2  =  A2(p,s,k)x2  +sB2(p)ze,  (34) 
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where 


and 


A2(fi,e,k)  = 


(  ea§2(ju)  +  ^2(e)E2  ea%z(n)  ae2i{p)^  2{e)\ 
-k2 

\MTddS2  Z2-£ae2M~2 


-k 

0 


0 

S2 


BoM  = 


(  jeM  \ 


c2=[l  0  0 


Proposition  IV. 1  There  exists  k*  >  0  such  that  for  any  given  k  >  k*  the  matrix  A2{(x,  e,  k )  is  Hurwitz  for 
all  n  G?  and  all  e  £  (0, 1] . 


The  proof  of  the  result  uses  elementary  arguments,  and  is  omitted  for  the  sake  of  brevity. 


Once  a  value  k  >  k*  has  been  fixed,  the  occurrence  of  k  from  the  arguments  of  A2(-)  is  removed  to  keep 
the  notation  streamlined,  and  to  stress  the  fact  that  the  selection  of  the  gain  k  has  already  been  made  and 
A2(n,e)  is  Hurwitz.  The  next  proposition  follows  directly  from  [19,  Lemma  2.1]. 

Proposition  IV. 2  There  exists  0  <  £*  <  1  such  that  for  all  e  £  (0,£*]  there  exists  a  solution  L(/x,  e)  of  the 
quadratic  matrix  equation 

B2(h)  +  A2(n,e)L (n,e)  -  eL(n,  e)Fe{p)  -  eL(p,,e)J12{p)C2L(h,e) 
which  is  approximated  as 

s)  =  -A21(h,e)B2 (h)  -  s  Af1(ia,£)B2(fi)[Fe(iu)  -  Jf2(n)C2AfL(n,£)B2(ii)]  +  0(e2)  . 

Note  that,  since  the  matrix  B2(/T)  appears  in  (34)  scaled  by  e,  the  change  of  coordinates  required  to  bring  (34) 

in  the  so-called  “standard  actuator  form”  is  given  by  X2  =  X2  —  eT(/x,  £)ze  ,  where  £  £  (0,  £*].  This  transfor¬ 

mation  yields  a  singularly  perturbed  system  in  the  form 

5e  =  Fe(n,£)ze  +  J1MC2X2 

£X2  =  A2(^,£)x  2  (35) 

where 

Fe(/x,£)  :=  Fe(p)  +£jf2([r)C2L(n,£) ,  A2(/i,e)  :=  A2(h,e)  -  e2 L(nt£)Jf2(n)C2  , 

from  which  it  is  now  obvious  to  conclude  the  following  result: 

Proposition  IV. 3  There  exists  £**  £  (0,£*]  such  that  the  system  (35)  is  asymptotically  stable  for  all  (/i,  e)  £ 
PX(0,£**]. 

The  main  result  of  the  section  is  then  summarized  as  follows: 

Proposition  IV. 4  Consider  the  dynamic  error-feedback  controller 
£2  =  pB  2f2  +  ^2fi 

u2  =  [Ai  +  4'2(p)£2  -  pk]£i  +  p^2(p)[kl  +  £2]6  -  pk2e2  .  (36) 

where  Ai  is  the  eigenvalue  of  the  dynamic  extension  (22),  k  >  0  and  p  >  0  are  tunable  gain  parameters,  and 
(£,E)  and  \k(p)  are  given  in  (28)  and  (30),  respectively.  Choose  k  >  k *  according  to  Proposition  IV. 1.  Then, 
there  exists  p*  >  0  such  that  for  all  p  >  p*  the  controller  (36)  applied  to  system  (25)  achieves  boundedness 
of  all  trajectories  and  asymptotic  regulation  of  e2{t). 
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3.  Outer-loop  Controller 


The  theoretical  development  of  the  outer-loop  controller  follows  along  similar  lines  as  the  inner  loop  and  will 
therefore  be  omitted  for  reasons  of  space,  and  only  the  final  form  of  the  controller  will  be  given. 

The  reference  model  for  the  velocity  reference  is  assumed  to  be  the  same  as  the  one  for  the  altitude 
reference  (obviously,  the  set  point  will  be  different).  Letting  Si  be  defined  in  the  same  way  as  <S2,  he. 
Si  =  S2,  an  internal  model  for  the  steady-state  control  for  the  velocity  dynamics  is  built  as  follows:  Let 
(Ei,  Hi)  be  the  pair  in  controller  form 


( 0 

1 

0  > 

(0\ 

Ei  = 

|  0 

0 

1 

5  Ei  = 

b 

\-0\3 

-0-4 

-05/ 

w 

where  the  coefficients  of  the  Hurwitz  polynomial  psi(s)  =  s3  +  05S2  +  oqs  +  < 7 3  are  chosen  arbitrarily,  in 
such  a  way  that  the  roots  of  PEi(s)  are  distinct  from  those  of  the  polynomial  ps1(s).  Let  Mi  be  the  unique 
nonsingular  solution  of  the  Sylvester  equation 


Mi  Si  —  E1M1  +  ziiQi  (38) 


and  let  T  :=  QiM1  1 .  A  candidate  internal  model  for  the  steady-state  input  u*(t)  is  constructed  as  the 
dynamical  system 


£3 

—  Ei£3  +  Hiei 

(39) 

Vi 

=  4>i^3 

while  the  overall  control  input  is  chosen  as 

Ui 

=  -fciei  +  Ti^3 

(40) 

where  ki  >  0  is  large  enough. 


E.  Controller  implementation 

The  combined  structure  of  the  controller  may  now  be  written  as 

£ 0  =  7^0  +  (P  +  i)Qf 

Cl  =  — ^Cl  +  Se 

1 2  =  ps  2C2  +  e2Ci 

1 3  =  E1C3  +  Si  (Vt  —  ri) 

Sc  =  -«Co  -nQf 

Se  =  [Ai  +  T2(p)H2  -  pk\ii  +  p'&2(p)[kl  +  E2]6  -  pk2(h  -  r2) 

(f)  =  -ki(Vt  -  n)  +  TiC3  ,  (41) 

whereas  the  procedure  to  tune  the  gains  of  (41)  is  summarized  as  follows: 

•  Choose  the  desired  damping  ratio  and  frequency  ((,w„)  for  the  reference  model. 

•  Choose  stable  polynomials  p^2  ( s )  =  s3  +  cr2s2  +  ays  +  00  and  (s)  =  s3  +  a5s2  +  oqs  +  <j3,  in  such 
a  way  that  the  roots  of  pz2  are  not  too  far  from  the  ju>- axis  (recall  that  the  roots  will  then  be  pushed 
to  the  left  by  increasing  p),  and  such  that  they  do  not  overlap  with  the  roots  of  Ps2(s)- 

•  Choose  the  position  of  the  rate  gyro  sensor  and  the  auxiliary  controller  parameters  7,  (3  and  k  such  that 
the  original  zero  dynamics  (18)  are  stabilized  for  a  given  range  of  fuel  conditions  and  flight  conditions. 

•  Choose  A  >  0. 

•  Fix  p  =  1,  and  compute  (the  smallest)  k *  such  that  the  matrix  A2(p.,  1,  k)  is  Hurwitz  for  all  p  considered 
(that  is,  for  all  possible  perturbed  models).  Fix  a  value  k  >  k*. 

•  Increase  p  until  the  system  (34)  is  asymptotically  stable  for  all  perturbed  models. 

•  Choose  ki  such  that  the  overall  closed-loop  system  is  asymptotically  stable. 
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V.  Simulation  Results  for  the  Robust  Output  Feedback  Design 


The  performance  of  the  robust  output  feedback  control  design  is  evaluated  on  the  nonlinear  model  (1). 
The  references  are  filtered  steps  with  respect  to  the  trim  condition  of  1000  ft/s  and  10000  ft  for  the  vehicle 
velocity  and  altitude,  respectively.  Since  the  airbreathing-hypersonic  vehicle  is  not  intended  for  highly 
aggressive  maneuvers,  we  shall  allow  ourselves  to  use  a  reference  that  settles  in  about  140sec  (1%  settling). 
Note  that  in  previous  works4-7  a  faster  reference  was  used,  settling  after  around  100  sec,  but  in  those  cases 
the  non-linear  model  did  not  encompass  heave  coupling  with  the  flexible  dynamics,  and  perfect  full  state 
measurement  was  assumed.  The  gains  for  the  robust  controller,  shown  in  Table  5,  were  tuned  and  validated 
by  verifying  that  they  stabilized  the  closed  loop  system  using  a  given  family  of  linearized  models  of  the 
nonlinear  plant.  Linearizations  have  been  performed  using  5  trim  conditions  and  two  fuel  conditions,  namely 
90%  and  50%  fuel  load.  Figures  13-14  show  the  results  of  closed-loop  simulation  for  a  90%  fuel  condition 
while  Figures  15-16  show  the  same  for  a  50%  fuel  condition.  It  is  evident  from  Figures  13(c)-13(b)  and 
15(c)-15(b)  that  the  tracking  performance  is  near  perfect  and  Qf  is  well  regulated  for  both  fuel  conditions. 
Furthermore,  Figures  13(a)  and  15(a)  show  that  the  inputs  are  fairly  well  behaved  for  both  fuel  conditions, 
as  they  are  within  their  saturation  limits  of  [0, 1]  for  0,  and  ±0.52rad  for  5e  and  Sc.  Here,  some  differences 
between  the  responses  for  the  two  fuel  conditions  can  be  observed,  as  a  smaller  </>  is  needed  to  maintain  the 
target  speed  at  the  50%  fuel  condition.  When  the  filtered  step  reference  is  applied,  a  noticeable  oscillating 
transient  due  to  the  coupling  with  the  flexible  dynamics  can  be  observed;  this  behavior  stems  from  the  use  of 
a  fairly  high-gain  feedback,  and  may  be  mitigated  by  a  more  careful  judicious  choice  of  the  controller  gains. 
The  angle  of  attack  deviates  only  slightly  from  trim  and  is  quite  well  behaved  for  both  fuel  conditions  as 
can  be  seen  from  Figures  14(a)  and  16(a).  Figures  14(b)  and  16(b)  show  that  the  flight  path  angle  is  kept 
reasonably  small  and  smooth  for  both  fuel  conditions.  Finally,  Figures  14(c),  16(c),  14(d),  and  16(d)  show 
that  the  flexible  states  and  controller  states  are  stable  and  settle  to  constant  values.  The  controller  dynamics 
were  implemented  in  a  balanced  realization  to  avoid  numerical  problems  with  the  control  canonical  form. 
Thus  the  states  shown  in  Figures  14(d)  and  16(d)  do  not  relate  directly  to  the  controller  dynamics  states 
given  in  previous  sections  but  through  a  linear  transformation  not  given  here  for  reasons  of  space. 

VI.  Conclusion 

In  this  study,  we  found  that  linear  observer-based  output  feedback  can  be  successfully  applied  to  a  non¬ 
linear  air-breathing  hypersonic  vehicle  model  for  a  very  limited  range  of  operating  conditions.  Modeling  of 
rate  and  acceleration  sensors  which  include  flexible  effects  provided  insight  into  the  effect  of  sensor  place¬ 
ment,  which  led  to  the  development  of  a  performance  index  focused  on  improving  system  observability.  The 
locations  of  high  observability,  near  the  nose  and  tail  of  the  aircraft,  agree  with  previous  work  in  sensor 
placement  for  flight  vehicles.  Since  the  observer  is  constructed  with  linearized  equations  of  motion,  the  ob¬ 
server  is  only  valid  for  an  operating  range  around  the  trim  condition,  and  offers  poor  robustness  with  respect 
to  parameter  variations.  This  limitation,  inevitable  with  a  linear  observer,  has  been  overcome  by  designing  a 
robust  dynamic  output-feedback  controller  which  does  not  rely  on  asymptotic  state  reconstruction.  Despite 
the  fact  that  the  robust  controller  constitutes  a  significant  improvement  over  the  observer-based  design,  there 
is  still  room  for  improvement.  The  fact  that  part  of  the  controller  is  designed  using  high-gain  feedback,  forces 
the  reference  trajectories  to  be  rather  slow,  to  avoid  inducing  considerable  oscillations  in  the  inputs  during 
the  transient.  A  more  sophisticated  design  based  on  adaptive  control  or  gain-scheduling  may  alleviate  the 
need  for  such  high-gain  feedback,  and  yield  improved  tracking  capabilities. 


Table  5.  Controller  Parameters  &  Sensor  Placement 


Rate  Gyro  Placement  (from  front) 

Aux.  Controller 

Other  Controller  Param. 

26  ft 

At  10 

(3  5 

7  1 

A  1 

k  138.6 

P  1 

k\  1 
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(a)  Inputs  to  the  plant  (b)  Measured  Pitch  Rate 


(c)  Output:  Velocity  &  Reference  (d)  Output:  Altitude  &  Reference 

Figure  13.  Plant  inputs  and  outputs  using  a  90%  fuel  condition 


(a)  Angle  of  Attack  (b)  Flight  Path  Angle 


(c)  Flexible  States  (d)  States  of  the  Balanced  Controller 


Figure  14.  Angle  of  attack,  FPA,  flexible  states,  and  state  of  the  controller  for  a  90%  fuel  condition 
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Figure  15.  Plant  inputs  and  outputs  using  a  50%  fuel  condition 


Time,  sec 


(c)  Flexible  States 


(d)  State  of  the  Balanced  Controller 


Figure  16.  Angle  of  attack,  FPA,  flexible  states,  and  state  of  the  dynamic  controller  for  a  50%  fuel  condition 
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